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Abstract. We here extend the well known Positive Real Lemma (also known 
as the Kalman-Yakubovich-Popov Lemma) to complex matrix- valued general- 
ized positive rational function, when non-minimal realizations are considered. 
All state space realizations are partitioned into subsets, each is identified with 
a set of matrices satisfying the same Lyapunov inclusion. Thus, each subset 
forms a convex invertible cone, cic in short, and is in fact is replica of all real- 
izations of positive functions of the same dimensions. We then exploit this re- 
sult to provide an easy construction procedure of all (not necessarily minimal) 
state space realizations of generalized positive functions. As a by-product, this 
approach enables us to characterize systems which can be brought, through 
static output feedback, to be generalized positive. 



1. Introduction 

For a half of a century, the Positive Real Lemma (also known as the Kalman- 
Yakubovich-Popov Lemma) has been recognized as a fundamental result in System 
Theory. We here extend and exploit it in various ways. Let C+ and C_ be the open 
right and left halves of the complex plane respectively, and P&, (P&) be the set 
of all k x k positive definite (semidefinite) matriceqj- Recall that a p x p- valued 
function F(s), analytic in C+ is said to be positive if 

(1.1) F(s) + F(s)* eP p sGC+. 

The study of rational positive functions, denoted by V, has been motivated from 
the 1920's by (lumped) electrical networks theory, see e.g. [7], [11]. From the 1960's 
positive functions also appeared in books on absolute stability theory, see e.g. [43] . 

m- 

A p x p- valued function of bounded type in C+ (i.e. a quotient of two functions 
analytic and bounded in C+) is called generalized positive QV if 

(1.2) F{iuj) + F(iu})* G P p a.e. wel, 
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where F(iuj) denotes the non-tangential limi10 of F at the point iui. 
Generalized positive functions were introduced in the context of the Positive Real 
Lemma (PRL), see |B] and references therein^. Applications of QV functions to 
electrical networks appeared in |35j . and to control in [40] . where they first casted 
in a Linear Matrix Inequality (LMI) framework, see e.g. [12) for more information 
on LMI. For more application of the generalized PRL, see [30] 

Both function sets V and QV are closed under positive scaling, sum and inversion 
(when the given function has a non-identically vanishing determinant). Yet these 
spaces have quite different properties, as we now illustrate. Let C pxp (X) be the 
space of C pxp -valued rational functions. Recall that a function "J" 6 C pxp (X) 
belongs to QV if and only if it can be factorized as 

(1.3) *(«) = G(s)P(s)G(-s*r, 

where G, P G C pxp (X), and P G V. Factorization of this nature appeared e.g. 
in [17] and [20], see also Observation 17. II below. The significance of (|1.3j) to scalar 
rational QV functions was recently treated in [4] and [5]. 

We now consider properties of the sum of two rational QV functions (series connec- 
tion in electrical engineering jargon). From (|1.2[) it follows that this sum is again 
in QV, but both the McMillan degree and the number of negative squares (roughly, 
the number of poles in C+) increase. For more details on the number of negative 
squares see [35] and [37]. Recall that a rational function "J is in QV if and only if 
the kernel 

*(s) + V(w)* 
s + w* 

has a finite number of negative squares in its domain of definition in C+. The 
number of negative squares of the sum of two elements in QV is preserved if, for 
instance, in JUJ *i(s) = G(s)P 1 (s)G(-s*)* and * 2 (s) = G(s)P 2 (s)G(-s*)* with 
Pi,P2 G V and the same function G G C pxp (X), see [5J Section 3] for the scalar 
case. 

In contrast, if one takes a state space realization sum, the McMillan degree of the 
resulting function does not increase, but it may turn to not generalized positive at 
all, see Example 17.21 below. One of the results of this paper is a partitioning of 
QV functions to subsets, denoted by QVir, v,p), closed under state space addition, 
while both the McMillan degree and the number of negative squares do not increase. 
See Theorem 17.31 below. 

We resort to some preliminaries. Let "J G C pxp (X) be of McMillan degree q, 
and analytic at infinity, i.e. lim &(s) exists. Namely, ^ admits a state space 

s — ► oo 

realization 

(1.4) *(s) = C(sI-A)~ 1 B + D L:= 

with A G C nxn , n > q, B,C* G C nxp and D e C pxp , namely, L G C^+rt x (n+p \ 
If the McMillan degree of ^(s) satisfies q = n, the realization is called minimal. 



•^This limit exists almost everywhere on iK because F is assumed of bounded type in C+ , see 
e.g. [22]. 

^The original formulation was real. The case we address is in fact generalized positive and 
complex, but we wish to adhere to the commonly used term: Positive Real Lemma. 
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We can now state the Positive Real Lemma (PRL) as presented in p~9l Theorem 1] 
(up to substituting the real setting by a complex one) 

Lemma 1.1. Let ^f(s) be a p x p-valued rational function in (|1.4[) and assume 
that q = n. 

(I) * G QV if and only if 

(1.5) HL + L* H = Q G P„ +p , 

for some H = diag{_ff, I p }, where H G <C nxn is Hermitian non-singular. 

(II) If * EV then in part (I) —H G P n • 

The aim of this work is to first extend this result to the non-minimal case, and then 
to use it to obtain a straightforward construction of all (not necessarily minimal) 
state space realization of V and QV rational functions. This is then used to 
describe the already mentioned partitioning to sets of the form QV(r, v,p) and to 
obtain other results, described below. 

The outline of the paper is as follows: The paper is composed of eight sections 
besides the introduction. In Section[5]we give a short review of the literature, which 
should be paralleled with a complementary survey we offered in our previous paper 
[4j. Our aim is not to provide a complete survey, but to raise, through samples, 
the intriguing observation that although the PRL has been a standard textbook 
material from the 1970's, see e.g. [7J, [UJ Chapter 3], 031 Section 4.4, Appendix] 
and [351 Section 8.5], it is not straightforward to cover the relevant literature. In 
simple words, there is too little of cross-referencing Q. 

In Section [3] an algebraic Riccati inclusion associated with necessity part of the 
generalized positive real Lemma is addressed. Sections |3J and [5J are devoted to 
showing that an algebraic Lyapunov inclusion associated with the sufficiency part 
of the generalized positive real Lemma is independent of the minimality of the 
realization. Some background material concerning sets of Lyapunov inclusions is 
reviewed in Section [6J In particular we provide a convenient parameterization of all 
matrices L satisfying the Lyapunov inclusion (jl.5|) where H is fixed and Q varies 
over all P. This Lyapunov inclusion formulation is then employed in Section [7J to 
provide a straightforward parameterization of all state space realization of QV 
rational functions. This allows us to describe QV functions as a union of replicas 
of positive functions. As an application of theses sets, in Section [8] we characterize 
all rational functions (vanishing at infinity) which can be rendered QV, through a 
static state feedback. Concluding remarks are given in Section [5] 

2. A HISTORICAL PERSPECTIVE 

We here review some of the relevant existing literature. 

As mentioned, the above version of Lemma 11.11 is from [19] and was repeated in 
[3D] . It was originally proved in [B]. A special case was later treated in [31]. The 
positive function case (part II) is well known and sometimes is referred to as the 
Kalman-Yakubovich-Popov Lemma and is dated to the 1960's. For an early full 
account see e.g. [7j Chapter 5]. An easy-to-read historical perspective is given in 
[121 Sections 1]. 



For example, it seems that 1191 was hardly ever cited. 
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A matrix formulation, see (|1.4|) . of the PRL was introduced in the PhD. thesis of 
P. Faurre, see [24j Theorem 4.2] and then in a book he co-authored, [25] Theorem 
3.1]. In fact it implicitly earlier appeared in [48] . The formulation through the 
Rosenbrock system matrix L (II. 4[) (for not necessarily positive systems) explicitly 
introduced in |19) and subsequently in [401 Lemma 8]. An interesting special case 
was studied in [27] Theorem 4] . The notion of Linear Matrix Inequality (LMI) was 
introduced in |48) . |40| was one of the early works recognizing the applicability 
of LMI framework to the (generalized) Positive Real Lemma (PRL), see also [T^J 
Section 2.7.2]. A comprehensive survey of the LMI approach to the PRL appeared 
in [29]. Unfortunately, in spite of its admirable reference list (201 items), some 
important relevant results are missing. 

Following (jl.ll) Positive functions map C+ to L(7), the set of matrices with a 
non- negative Hermitian parlQ. Analogously, following (jl.2l) . a QV function maps 
iR (after removing all poles of the function) into L(7), see Observation 17. II below. 
Closely related function sets are addressed in the literature: 

• Bounded functions mapping C + to weak contractions while generalized 
bounded functions map tSL to weak contractions. Versions of the PRL for 
bounded functions appeared [23] Theorem 4.2] [7] Section 7.2] and for gen- 
eralized bounded in Q32 Equation (6)]. An interesting subclass is treated in 
[TJ Theorem 2.1]. An important result on a subclass of generalized bounded 
functions where iR is substituted by R, appeared in [28] Theorem 3.2]. 
See also [2] Theorem 2.12] and [38] Chapter 21]. 

• Caratheodory functions, analytically mapping the open unit disk to L(/) 
and generalized (=pseudo) Caratheodory functions, mapping the unit circle 
(from the poles of the given function have been removed) to L(/). Versions 
of the PRL for Caratheodory functions appeared in [3S] . 

• Schur functions, analytically mapping the open unit disk to weak contrac- 
tions and generalized Schur functions, mapping the unit circle to weak con- 
tractions. Version of the PRL for Schur functions appeared in [3] Theorem 
2.5], [47] Section II], [49] and for operator valued in [10]. For generalized 
Schur functions see [TH1 Theorem 2] and for an interesting special case see 
[27l Theorem 3] . A time- varying extension of the PRL to generalized Schur 
functions along with a thorough study of various applications is are given 
in [SHI Theorem 1.2.5 and Appendix 3A]. 

• Nevanlinna functions, analytically mapping the open upper half plane to 
iL,(I) and the generalized Nevanlinna functions map M. to see 
e.g. [17] and |20) . Minimal realization of infinite dimensional generalized 
Nevanlinna functions was studied in [21] . 

As already stated, we do not aspire to provide a survey of PRL related results, and 
we are aware of additional references dealing with the subject, not mentioned here. 
We focused on the scattered nature of the literature related to the generalized case. 

3. Generalized positive lemma necessity and the Riccati equation 

It has been long recognized that with the part (b) of Lemma 11.11 (dealing with 
positive functions) one can associate an algebraic Riccati equation, see e.g. [7J 

^The notation L(7), with L honoring A.M. Lyapunov, will be formally defined in I l3.ll ! below. 
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Section 5.4] [25j Chapter 5] and [12| Section 2.7.2]. For Schur functions version 
see e.g. p~8| Theorem 2.1]. We now address the extension of this result to QV 
functions. It first appeared in the context of a variant of Generalized Bounded 
functions (where iM. was substituted by K) in [28} Theorem 3.2] and in [2j Theorem 
2.12]. For another variant of this result, see [38l Section 20.1]. In Proposition 13.21 
below we provide a simple proof of the result using the system matrix formulation 
(|1.4p . employed all along this work. We shall find it convenient to resort to the 
following notation of sets of matrices sharing a common Lyapunov factor: For a 
r x r Hermitian non-singular matrix H, define the sets of r x r matrices, 1<(H) 
and L(-ff) as, 
(3.1) 

(a) TL(H) := {L : HL + L*H G P r } (b) L(H) := {L : HL + L*H G P r }. 

In particular, (L(7) ) L(J) is the set of matrices with positive (semi) definite 
Hermitian part. 

Proposition 3.1. Consider (|j~4) . (TTT5]) with H = diagjtf, I p } and H G C nxn 
Hermitian. Assume in addition that D G L(/ p ). Let us define the following n x n 
Riccati expression, 

M :=H{A-B(D + D*)- 1 C) + (A — B(D + D*)~ l C)*H 
( ' ' -HB(D + D*)- 1 B*H - C*{D + D*)~ l C. 



Then Q G P n + P , if and only if in 
(3.3) M G P„ . 

Proof Using the fact that D G L(/ p ) one can employ the classical Schur's 
complement, e.g. 33, Theorem 7.7.6], to write down Q in (jl.5l) explicitly, 



Q 



HA + A*H HB + C*\_fl n R\(M \ f I n 
C + B*H D + D* ) \0 I p ) V D + D*) \R* I p 



where R := (HB+C*)(D+D*)~ 1 and M is given in (J3HJ). Thus indeed Q G P n+P 
if and only if M G P„. □ 

We can now re- formulate the necessity part of Lemma 11.11 

Proposition 3.2. Let ^ G QV be a p x p-valued rational function so that 
lim ^(s) := f exists. Furthermore assume that D G L(/ p ). Let q be the 

s — > oo 

McMillan degree of ^. 

(TJ Then, ^(s) admits a state space realization (11.4[) . so i/iai if is minimal (q = n) 
and it satisfies the Riccati inclusion in (|3.2j) . p.3j) /or some nx n Hermitian non- 
singular H . 

(II) If * £P t/ien m part (7j -F G P„. 

Note that the technical condition D G L(/ p ) in Proposition ^ . 2l is indeed restrictive. 
For example, many system of interest have a zero at infinity and thus are excluded 
from the discussion. On the other hand, whenever Vl'(s) = C(sL — A)~ 1 B + D is a 
QV function, from the above discussion it follows that D G L(J p ), hence one can 
always construct another QV function ^>(s) = C(sL — A)~ 1 B + D, with D G L(J p ) 
so that e > — D\\, where e > is arbitrary. 
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4. Generalized positive lemma sufficiency - an extension 

The sufficiency statement of the (generalized) Positive Real Lemma (PRL) of matrix 
valued rational functions was first proved in [B] , under the assumption of minimality 
of the realization {q = n). We now address the question of relaxing this minimality 
constraint. This problem was treated in the framework of positive functions in 
[49l Lemma 6] and in the framework of functions satisfying D £ ^(I P ) (as in 
the previous section) in [38l Section 21.3]. In a different formulation see also [46j 
Theorem 1]. In [55] a proof of the sufficiency statement of the (generalized) PRL, 
removing the minimality of realization condition, was presented. Unfortunately a 
(redundant) spectral condition on A was imposed there. The result of Proposition 
14.21 below, avoids any restriction. 

In addition, in Proposition 14.21 below we show that one can bound the number of 
poles of ^(s) in each open half plan^. To this end, we need some preliminaries. Re- 
call that for a matrix A £ C" x ™ one can associate a triple: inertia(j4) = (u, 6, n), 
with v + 5 + n = n, if A has v eigenvalues in C_, tt eigenvalues in C+ and 
S eigenvalues on iR, see e.g. pSI 2.1.1]. Let A, H £ C nxn with H Hermitian, 
be with inertia, 

(4.1) mertia(A) = (ua, 6a, ka) mertia(.ff) = (u, 0, n — v) v £ [0, n], 
i.e. H is non-singular. Consider now the Lyapunov equation 

(4.2) HA + A*H = Q eP„ . 

Recall that from a pair A £ C nxn , B £ C nxp , the following controllability matrix 

may be constructed C := [B : AB I • • • : A n - 1 B}. Then X cont (A,B), the control- 
lable subspace associated with the pair A, B is given by the range of C and 
-^cont(A B) , its orthogonal complement, is given by the null-space of C*, see 
e.g. [31 Definition 2.4.8]. Similarly with a pair A £ C nXn , C £ C pxn , one can 
associate a observable subspace X b s (A, C), given by X b s (A, C) — X cont (A* ,C*) 
(and Xob^AC)- 1 =A ront (A*,C*)- L ). For A, Q in @3) denote, 

(4.3) m := dim X ohs (A, Q) x . 

Namely, pair A, Q is observable whenever m = 0. We can now cite the following 
important result of R. Loewy [39] , adapted to our framework, 

Theorem 4.1. Let A,H and m be as in (|4T|) . (j4~2"]) and (|Pjt . Then, 

v > va ^ max(0, v — m) n — v > tta > max(0, n — v — m). 

We can now state the extended sufficiency part of the PRL. 

Proposition 4.2. Let "J in (|1.4j) be a p x p-valued rational function of McMillan 
degree q, n > q. Assume that state space realization in (11.41) satisfies the Lyapunov 
equation (|1.5[) with H = diagjiJ, I p }, where H is n x n Hermitian. 

(I) If inertia(-ff) = (y, 0, n — v), for some v £ [0, n], then, ^ is a QV function 
with at most v poles in C_ and n — v poles in C+ . 

(II) If in part (I) -H £ P„, i.e. v = n, then 



Although in a different framework, bounds of a similar nature can be found in |28l Theorem 
3.4] and subsequently in [2] Theorem 2.12] and in [381 Section 21.2]. 
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Proof : I. Indeed assume that (|1.5p is satisfied with H = di&g{H,I p }, H 
Hermitian nonsingular and L as in (|1.4j) . Note that in ()1.5|) Q is in P„+ p , thus 
its upper left block is in P„. Namely, (|4.2|) is satisfied, so by Theorem 14. II the matrix 
A has at most v eigenvalues in C_ and n — v eigenvalues in C+. Recall that 
the poles of ^ are determined by the eigenvalues of A, see (|1.4p . 
Next denoting S :— diag{— sl n , P } we note that 

(4.4) HS + S*H G Pn+ P , 

for all s G iR. Take nowQ L:=L + S = (' sI " c + A ■ Then for all s G iM. also, 

(4.5) ffX + L*H 6 P Il+p . 

Next, recall that for arbitrary constant matrix ^ G QV, if and only if (^ + T) G ^T 5 , 
for arbitrary constant matrix — T* = T G C pxp . Thus, up to a shift by a skew- 
Hermitian matrix, we can assume that vE'(s) in (ll.4[) is almost everywhere invert- 
ible. Thus, a straightforward calculation (see e.g. [33l 0.7.3]) results in, 

\sl - A)- 1 (B*(s)- 1 C(s/ - A)- 1 - I) (si - A^B^is)- 1 " 
*(s)- 1 C(s/- A)' 1 f(s)- 1 

Multiplying (|4.5[) by (i*) from the left and L^ 1 from the right, yields 



(4.6) i/i- 1 + (L~ L )*H e ¥ n+p 

for all s G «R. Now in particular the p x p lower right block of (|4.6[) satisfies, 

(4.7) ^(s^ + r^s)- 1 )* = (Hi- 1 + {L- 1 )*H) GP P 



' 22 

for all s G iR. Thus, 4'(s) _1 is in and hence also ^(s). Thus, the first part 
of the claim is established. 

(b) Positive functions. If —H G P„ the relation in (|4.4[) holds for all s G C+ and 
subsequently, also (|4.6[) and (|4.7|) . Hence, $ e P, so the proof is complete. □ 



In the next section we scrutinize some aspects of Proposition [ 

5. NON-MINIMAL REALIZATION AND BOUNDS ON INERTIA - A CLOSER LOOK 

Roughly, application of Theorem 14. II to Proposition 03] suggests that the "further" 
from minimality the realization is, the cruder is the bound on the number of poles 
in C+. We here illustrate the "at most" clause in the statement Proposition 14.21 
with respect to the number of poles in each open half plane. 

Example 5.1. All functions considered in this example are so that in (|1.5[) H = diag{-ff , 1} 
with H = diag{ — 1, 1}, see also (|4.2[) . Namely the corresponding rational functions 
have at most one pole in each open half plane. We present the rational function 
along with the corresponding system matrix. 

/0 1 1\ 1-1 0\ l-\ 1 

L a = i i o La = oil L 7 = o i o 

\1 0/ \ 1 1/ \ 1 

f-l -1 l\ /0 l\ /-l -1 2\ l-\ 1 -2 

Lg=[ i i i Le = o i o L„ = -i i i £0=1 1 1 

1 -1 0/ \1 1/ V 2 1 1/ V-2 1 1 



'It is interesting to note that the zeroes of *S>(s) are the points s for which L is singular. 
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1>s(») = w M s ) = M s ) = Ms) = ■ 

Indeed, ip a (s) and ip v (s) = ipg(s) have a pole in each open half plane. ipff(s) has 
a single pole in C+, while ip^ (s) a pole in C_. The poles of ips( s ) an d ^(s) 
are confined to the imaginary axis. □ 

We now point out that unlike to the approach of [23 Theorem 2], in the PRL 
framework, see Proposition 14.21 there is no restriction on the spectrum of L in 
(|1.5p . nor on the spectrum of its upper left block, A in (|4.2[) . Formally, this 
follows from Theorem 14.11 Intuitively, restrictions of the form Xj + A£ ^ are 
necessary when Q, the right hand side of the Lyapunov equation is given and 
a unique Hermitian solution is sought, see e.g. [3H Corollary 4.4.7]. In contrast, 
here we address a genuine Lyapunov inclusion. In fact, in each of the cases Lp, 
L 7 , Lg, L^, L^ and Lg in Example 15.11 the A matrix (the n x n upper left 
block) does not satisfy this spectral condition. 

6. Convex invertible cones and the Lyapunov inclusion 

In this section we explore some properties of the set L(ff ) (|3.1[) . to be used in the 
sequel in conjunction of the PRL, see e.g. Theorem 17.31 

We next resort to some background. Recall that a set of a square matrices is 
called Convex Invertible Cone, cic in short, if it is closed under positive scaling, 
summation and inversion. More precisely, it may include singular elements provided 
that the inverse of every non-singular element, belongs to the set. Thus, the set 
P is a cic. For a more detailed study of cics see e.g. [T3J Section 2], [14l Section 
2] and Sections 2, 3] Recall that for an Hermitian non-singular matrix H we 
defined in (|3.ip sets of matrices sharing a common Lyapunov factor, h(H), h(H). 
The following properties of these sets are fundamental to our discussion. 

Theorem 6.1. Let H be rxr Hermitian nonsingular, i.e. inertia(if) = (v, 0, r — v) 
for some v G [0, r\. 

(i) L(if) is an invertible cone of matrices sharing the same inertia as H . It 
is a maximal open convex set of nonsingular matrices^. 

(ii) L(_ff) is a closed invertible cone of matrices with at most v eigenvalues 
in C_ and at most r — v eigenvalues in C+. It is a maximal convex set 
with this property. 

(iii) 11(H)* = £(#-_]■). 

In particular. L(_ff)* = L(-ff), if and only if, up to positive scaling, H is 
an involution (H 2 = I) . 

(iv) Let E be an involution which commutes with H then, 
Eh(H) = L(H)E = h(EH). 

(v) For arbitrary Hermitian involution E, Elu(E) = li(E)E = Ij(7). 

Proof (i) See [131 Proposition 3.7], [T51 Observation 2.3.3]. 

(ii) Part of the claim appeared in [14j Proposition 2.4(1)]. To establish the inertia 
property assume first that H = I r . By Theorem 14.11 all matrices in L(J) have 
inertia (0, r — it, it) for some tt G [0, r\. Let B be an rxr matrix 
not in L(J), namely the Hermitian part of B has a negative eigenvalue, i.e. 
min \ k (B + B*) = -(3 for some f3 > 0. Take now A = §1 + \(B* - B). 



^For an impressive particular converse, see [8], [9]. 
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It is straightforward to verify that A E L(7), but A + B = § I + \{B + B*). 
Thus, min \k{A + B) — —4, i.e. in contrast to all elements of L(I), the 

k— 1, ... , r 

matrix A + B has eigenvalues in C_ so this part of the claim is established for 
H = I. To address a general Hermitian non-singular matrix H , exploit the fact, 
see Observation 16.31 below, that one can always find a non-singular V and an 
involution E v §3J so that E v V- l Z(H)V = V~ 1< L(H)VE V = L(I), so the claim 
is established. 

(iii) The claim for arbitrary H appeared in [T31 Equation (3.6)]. If H is a scaled 
involution, i.e. H^ 1 = aH for some a > 0, the claim is follows from multiplying 
h(H) in dm]) by H^ 1 from both sides. 

For the other direction assume that (L(H)) = h(H). This relation is invariant 
under unitary similarity. Thus, without loss of generality assume that H is (real 
non-zero) diagonal, say H = diag{/ii , . . . , h r }. Take now A to be equal to H, 
except a single non-zero off diagonal element x at the location jk, where j > k, 
i.e. A is lower triangular. A straightforward calculation shows that A E h(H) 
is equivalent to 2\hj\ > \x\. By assumption, also A* E L(i?), which in turn is 
equivalent to 2\hk\ > Thus, hj — ±hk and since true for all j, k this claim is 
established. 

(iv) This claim is proved for h(H), in [T3J Lemma 3.6], see also [151 Proposition 
3.2.2]. The case h(H) is similar and thus omitted. Item (v) is in fact a special 
case of item (iv). □ 

Note that item (iv) in Theorem 16.11 says that if L E h(H ) , for some Hermitian 
non-singular H, then both matrices EL and LE belong to h(EH), whenever 
E is an involution which commutes with H. Using this, along with item (v) in 
Theorem 16. 1[ we state the following. 

Corollary 6.2. For natural n,p and vE [l,n] let us denote H = diag{I l ,, —I n - V , I p }, 
Hi = diag{-/„, /„_„, I p } and H 2 — diag{-/„, I p }. Then, 

HoL{H x ) = U,H X )H = L(JT a ), 
H h(H 2 ) = h{H 2 )H = L(Hx). 

In addition, HjT,(Hj) = L(Hj)Hj = L(I n+p ) j = 1, 2. 

Namely, there is one-to-one correspondence between the sets L(_ffi), 1 L{H 2 ) and 
L(/ ?i +p). 

We shall find it convenient to introduce the following notation for I x I signature 
matrices, 

(6.1) E v j := diag'l-7, , !,_„} v E [0, I}. 

Whenever the dimension / is evident from the context we shall simply write E v . 
We can now cite the following known facts see e.g. [T31 Lemma 3.4]. 

Observation 6.3. Consider the r x r nonsingular matrices V and H = H* . 
Then the following relations hold in (|3.1j) . 

(a) V- l l.{H)V = h(V*HV) (6) V^L^V = Z(V*HV). 
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In particular, if inertia(-ff) = (y, 0, n — v), for some v G [0, r], V may be 
chosen so that, 

(a) V~ 1 l.(H)V = L(£„) (b) y- x L(if)y = T,(E V ), 

where E v is as in (16.11). 



From Observation [673] it follows that, up to similarity, (J L(_E„ jr ) covers all rxr 

!/=0 

matrices. We next point out that technically this is not a proper partitioning. A 
straightforward substitution in (|3.ip with both H — E v and E u+rj , reveals that 
these sets are "nearly" distinct. 

Corollary 6.4. Let v > and -q > 1 be so that r > u+r] and E v is as in (|6.ip . 

If L G {h(E„) P| L(£' ;y+ „)} </ien L = ( Ql o t 3 o | where Ti,T%, T3 are 

\ k* Q 2 + T 2 J 

skew-Hermitian matrices of dimensions v x v, r\ x ry and (r — v — rj) x (r — v — if), 
respectively, K G C I/x ( r - I/ -'j) j s arbitrary and the other blocks are so that the 
(r — 77) x (r — 77) matrix r 2 ^, 1 2 g 2 ^ is positive semi-definite. 



Observation 16.31 also suggests that for every non-singular Hermitian H the set 
L(7J) is isomorphic to L(7). We conclude this section by showing that in turn, 
one can describe L(J) as a sum of two cics: P and T, the set of skew-Hermitian 
matrices, see e.g. [151 Proposition 3.2.5(h)]. Furthermore, each may be described 
by the convex hull of its extreme directions. To this end, we need to introduce OP, 
the set of rank one orthogonal projections, i.e. 

OP = {-k = xx* : x G C r x*x = 1 }. 

Observation 6.5. /. Let E be as in (|6.1I) . then 

EL(I) = h(I)E = L(E), 
EL(I) = h(I)E = L(E). 

II. [HI Proposition 3.2.5] 

L(7)=P + T E(J)=P + T. 

III. The sets P and T may be constructed from orthogonal projection. Indeed, 

3"! : a 3 > 

where tti, ... , 7r r are all distinct. Similarly, 
where tt\, ... , 7r r are all distinct. 

To summarize, for arbitrary E there is a one-to-one correspondence between 
the sets h(E) and L(7). Thus, it suffices to construct the latter set. Indeed, 
L(7) = P + T. Now, P G P can always be parameterized by non-negative scalars 
ai, ... , a r and r distinct points on the || H2 unit sphere. Similarly, TgT 
can always be parameterized by real scalars pi, ... , p r and r distinct points 
on the || H2 unit sphere. 
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In fact, parameterization of a point on || ||2 unit sphere can be further simpli- 
fied. Note that it E OP may be identified with a point on the || ||a unit sphere 
{xeC : x*x = 1 }. Now, through polar coordinates there is a one-to-one corre- 
spondence between this unit sphere and { y 6 R r (''~ 1 ) ; 2n > ||y||oo}- For example 

/cosffli) cos(e 2 )e i,71 \ 

for r — 3, x = sinfex) cos(9 a )e i "= Qj,T1k S [0, 2%), j — 1,2 and fc = 1,2,3. Thus 

\ sm(0 2 )e i '>3 / 

elements in OP can be parameterized by points in the real "box" [0, 27r) r(r ~ 1 ). 

In the next section we identify the matrix cic h(H) with the set of system matrices 
associated with a subset of QV functions, see Theorem 17.31 below. These (not 
necessarily minimal) realizations, cover all QV functions with no pole at infinity. 

7. CONVEX INVERTIBLE CONES OF REALIZATIONS OF GENERALIZED POSITIVE 
FUNCTIONS OF PRESCRIBED PARAMETERS 

We start with a couple of known related results. First, recall that in (ll.2[) we have 
already described QV functions as map from iM. to L(7). We now formalize this 
fact. 

Observation 7.1. One can view p x p-valued QV functions as a cic of rational 
functions (almost everywhere) analytically mapping iR to h(I p ), (|3.1j) . and V 
as a subcic mapping C+ to L(ip). 

In [HI Section 4.4] analogies were drawn between the set V of scalar rational 
functions and the matrix set h(H ) for H 6 P. We now use Observation 17.11 to 
introduce an analogy between the sets QV and L(7). Recall that 

vl,{i)v* c E(J) 

(and if in addition V is non-singular, then, 

FQVF* C QV 

Recall also that 

As another association between the matrix cic L(JJ) and a cic of rational func- 
tions, we cite the following. Consider the set of all rational functions of the form 
C(sl — A)~ l B admitting balanced realization, i.e. 

HA* + AH = BB* HA + A*H = C*C, 

for some non-singular Hermitian H . Each of these sets forms a cic of state space 
realizations. In particular, this allows for simultaneous model order reduction of 
uncertain systems, see [UJ Section 5] for details. 

We now identify the matrix cic h(H), (|3.ip . with a cic of system matrices, (|1.4[) . 
associated with a subset of QV functions. As a motivation recall that the set 
of p x p QV functions is a cic of rational functions. However, the McMillan 
degree of a sum, is roughly the sum of the McMillan degree of the original functions. 
Now, if one considers a pair of p x p QV functions, of McMillan degree at most 
n, admitting state space realizations of the form (jl.4p , the sum of the respective 
system matrices is associated with a p x p rational function of McMillan degree 
at most n. However, this "sum" function may be not generalized positive. This is 
illustrated next. 



V e C" x ™, 
Vh(I)V* C L(J)). Similarly, 

FeJ. 
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Example 7.2. For simplicity, consider two (minimal) state space realizations of 
the scalar function ip e (s) = ips{s) — from Example 15.11 L$ =( i 



' 1 

1 -10/ 



and L e = ( i -i -i | (note that L t = V 1 L$V with V — diagjU , 1} where 



-i o 



V = (° y 1 ) ). Let now be a convex combination of these two realizations, 



i.e. — ^(Ls + L e ) =[i o oj. This is a (minimal) state space realization of 



o-io 



i})q{s) = ^rpi | = — lj which is not generalized positive. □ 
Coordinates transformation 

Following Proposition 14.21 we have focused our attention on a special case of the 
set L(iJ), where r = n + p and the Lyapunov factor H is block diagonal of the 
form H = diag{H, I p }, where H is n x n Hermitian non-singular. 
From Observation 16.31 it follows that up to similarity over "L(H)\ H _ d . {ft : } i one 

can confine the discussion to the case where in addition in (|3.ip H = E un , see 

(E3). 

Recall that coordinates transformation means that whenever V = diag{y,/ p } 
with V n x n nonsingular, with L in (II .4[) . V^LV is another state space 
realizations of the same rational function ^(s). Thus, taking in (|6.1[) I = r = n+p, 
without loss of generality we can focus on sets of (n+p) x (n + p) matrices of the 
form 

L(H) H = diag{£^„, I p } v e [0, n], 

partitioned as . We now find it convenient to denote by 

(7.1) QV(r, v, p) r>2 p e [1, r - 1] v&[0,r-p] 

the set of all p x p-valued rational functions obtained by (|1.4p and (|1.5p (recall, 
A is (7' — p) x (r — p) ). From item (ii) of Theorem 16.11 we have the following: 

Theorem 7.3. Given a set QV(r, is, p) as described in (17. ip . The set of all the 

corresponding state space matrices forms a closed invertible cone of realizations of 
functions with at most v poles in C_ and at most r — p — v poles in C+. 
Moreover, this is a maximal convex realization set with this property. 

Maximality is in the sense described in proof of item (ii) of Theorem 16.11 Recall 
that each of the realization matrices L associated with functions in QV(r, v, p) 
has at most v eigenvalues in C_ and at most r — v eigenvalues in C+. Now if 
L is a rxr realization matrix associated of a rational function ip$.GV(r, v, p), 
then one can always find L so that L + L has v + 1 in C_ or r — v + 1 
eigenvalues in C+. 

Before proceeding, we now find it convenient to denote by QV-mi^r, v, p) the 
subset of QV(r, v, p) where in addition the realization is minimal (i.e. the 
McMillan degree q is equal to r — p). Under this terminology the necessity 
part of the PRL, i.e. Lemma Tl. 11 and Proposition 13.21 is restricted to elements in 
GV min {r, v, p). 

Example 7.4. We next illustrate some aspects of the correspondence between the 
matrix set h(E v ) and the set of rational functions QV(r, v, p), along with its 
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subset GVraia{r> v, p). As before, we concentrate on the case where r — 3, v = 1 
and p = 1. 

(i) To illustrate Theorem 17.31 the functions ip a (s), tpp(s), ip-y(s), i^s{s), 
ip^(s) and ip v (s) from Example |5~T1 are all in GV(3,1,1). 

(ii) Starting from L, one may obtain, minimal and non-minimal realizations. 
4>p(s), ip.y{s) and ^{(s): from Example l5.1l are in QV(3 7 1, 1) \ GV m in{3, 1, !)• 
In contrast, starting from rational functions, ip^(s) and V'fl 5 ) can be writ- 
ten as part of GV mi n(2, 1, 1), i.e. L 7 = °) and i" 1 = L € = (° J) 
are the corresponding system matrices. This now conforms with Proposi- 
tion [776] part (I) below noting that they are positive. 

(iii) L(i?i) is a convex invertible cone, cic. Samples of it are all realization 
presented in Example 15.11 In particular, = L^ 1 and Lp = \(L^ + Lg). 
In contrast, in Example [7721 L e is constructed from Lg so that they do 
not share a common Lyapunov factor, of the form diag{/f , 1}, to the LMI 
in (fT31> . 

(iv) Even in 5'P m in(3, 1, 1) elements may not have a pole in each open half 
plane, see e.g. ips(s) in Example 15. II 

(v) In contrast to the set QV(r, v, p) , the set of realization matrices associated 
with the subset GV m i n (r, v, p) , is not convex. Recall that on the one hand, 
ipr)(s) = ipe(s) is in (3, 1, 1), i.e. L v , L e are the respective minimal 
realizations. However, Lp = \{L, n + Lg) is a non-minimal realization of 

□ 

We now use the Lyapunov inclusion formulation to establish a strong link between 
positive and generalized positive rational functions. 

Lemma 7.5. Let r, v, p with r > 2, p G [1, r — 1] and v € [0, r — p] be given 
and let J = diagj/^, — 7 r _ p _„} (i.e. J — —E v ^ r ^ v ). 

Let (t, be a system matrix associated a function in QV(r, u, p) and let 

be associated with a p x p-valued positive rational function whose realization is of 

dimension r — p. 

Then, the following is true. 

(-0 {^c J d) an d (cj d) are t wo T — p dimensional realizations of the same 

p x p-valued positive rational function. 
(I) \£ J ~ 3 ") an d a) are two realizations of the same function in QV{r, v, p). 

Proof : To see that these are two realization of the same rational function, recall 
that for arbitrary involution J, one has that, 

C(sl n - JA)- X JB + D = C(sJ - A)~ 1 B + D = CJ(sI n - AJ^B + D. 
For the equivalence, the claim follows from Proposition 14.21 and Corollary 16.21 with 
n = r — p and H Q = J. □ 

This can be formalized in the following stating that for given r > p > 1, the set of 
p x p-valued rational generalized positive functions whose realization is of dimension 
r — p, can be described as r + 1 — p replicas of its subset of positive functions. 

Proposition 7.6. Let r,v,p be arbitrary. 
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(I) If r = v + p and v > 1 then, 

gV(r, v, p) Ut CP. 

(II) The set QV(r, v, p) is state-space equivalent to the set of positive functions 
QV{r, r - p, p). 

(III) For given r > p > 1, the set of p x p-valued rational generalized positive 
functions whose realization is of dimension r—p, can be described as 

r—p 

U SV(r, v, p). 

To see that the inclusion in item (I) of Proposition 17.61 is strict, take for example 
ijjj(s) from Example 15.11 where r — 3, v = 1 and p = 1, i.e. it is a positive 
function in QV (3, 1, 1) \ ^^111111(3, 1, 1), but neither belongs to QV(3, 1, 2) nor 
to QV{3, 2, 1). See also item (ii) in Example El 

Strictly speaking the sets QV(r, v, p) do not offer a partitioning of generalized 
positive function to distinct subsets. Namely, from Corollary 16.41 is follows that 

QV(r, v\, p) V\QV(r, V2, p) ^ 0- For example, the system matrix L =(0 a f>) 

\0 -6 Oj 

with a > 0, is a realization of tp = -3- . As L e L(I) D h(Ei), this V belongs to 

gv(3,o,i)ngv(3,i,i). 

Observation 16.51 offered an easy-to-compue construction of L(J) and subsequently 
of L(£'), where E — E u>r and v € [0,r — p] is arbitrary. Theorem 17.31 and 
Eq. (|7.1|) identify the matrix set L(i?) with QV(r, v, p). Thus, through item 
(III) of Proposition 17.61 we here offer a construction of the set of all p x p-valued 
rational QV functions whose realization is of dimension r—p, where r > p > 1, 
are arbitrary. 

LMI control theory and Matlab based LMI procedures, were developed essentially 
for V functions, see e.g. [12] . However, as stated in Proposition 17.61 from real- 
ization point of view, for arbitrary r > 2, p G [1, r — 1] and v G [0, p — r\, the set 
QV(r, v, p) functions is equivalent to the set QV(r, r—p, p) positive functions. 
It is thus of interest to try to adapt some of the LMI theory to QV functions. 

To further motivate the introduction of the set QV(r, v,p), in the next section we 
provide a classical control interpretation of it. 

8. TURNING A FUNCTION GENERALIZED POSITIVE THROUGH STATIC OUTPUT 

FEEDBACK 

Let F(s) be a p x p-valued rational function vanishing at infinity ( lim F(s) = 0). 

s— f 00 

Thus, it admits a space realization of the form F(s) = C(sl — A)^ 1 B, 

(8.1) x = Ax + Bu, y = Cx, 

where u,y are p— dimensional and x is n— dimensional . Using (|1.4[) . the 
corresponding system matrix is L = {f, q\ Assume that a static output feedback 
is applied i.e. 

(8.2) u = Ky + u', 



9 The case D ^ is more involved and thus omitted. 
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where K is a constant p x p matrix and vl is an auxiliary input. Thus, the 
resulting closed loop system (mapping v! to y) is 

(8.3) F cl (s) := C(sl- A d )~ l B, A cl :=A + BKC. 

The associated system matrix is, L c \ = { A q $y Robust stabilization of a system 
of the form (|8.1I) by static output feedback (|8.2p is quite known, see e.g. [TBI 
Theorem 3.2], [23] Section 3] [31] [32] and [44]. It is well known that positive 
functions are closely related with robust stability, see e.g. [43], [24]. We here 
introduce a characterization of all systems which may be turned generalized positive 
through static output feedback. 

Proposition 8.1. Let F(s) be a p x p-valued rational function vanishing at 
infinity ( lim F(s) = 0) with a state space realization as in (|8.1j) . 

s— »oo 

(I) There exists a static output feedback, of the form of (|8.2j) . rendering the 
closed loop system, F c \(s) in (|8.3p generalized positive, if and only if, there 
exists a nonsingular Hermitian H G c™ x ™ so that the open loop system 
(|8.1[) satisfies, 

(a) C = -B*H, 

(b) v*(AH~ 1 -\- H~ l A*)v > for all vector v in the null-space of B* . 
(II) Let r,v,p, where v G [0, r—p] and p € [1, r — p] be given. There exists a 

static output feedback, of the form of (|8.2[) . so that F c \(s) in (|8.3p belongs 
QV(r,v,p), if and only if, up to coordinate transformation, the open loop 
system (18. 1|) satisfies, 

(a) C = -B*E V , 

(b) v*(AE v + E l ,A*)v > for all vector v in the null-space of B* . 

(III) Let r,v,p, where v G [0, r — p] and p G [1, r — i'] 6e given. The 
set QV{r,v,p) is invariant under static output feedback (|8. 2[) whenever 

Keh(-I p ). 

(IV) 7/m pari (I) —H G P n or in parts (II), (III), v = r—p, then the resulting 
closed loop function F c \ is in V . 

Proof: (I) Following Proposition S2] denote H = diag{#, I p } with H G C™ xn 
Hermitian nonsingular, W :— HL C \ + L* cl H and recall that having F c i G QV is 
equivalent to We P r 

UT—{HA + A*H C* + HB\, I HBKC + (SA'C)* H 0\ 
\ C + B*H P \ Op)' 

Thus, having W 6 P r , implies condition (a) in the claim. Substituting back one 
has that 

\\r = (HA + A'H 0\i I -HB(K + K*)B* H 0\ —(HWH 0'\ 
\ o o p y \ P ) \ Op) 

where W = AH^ 1 + h(- 1 A* - B(K + K*)B*. Now, W G V r is equivalent to 
W G P„ (recall r = n + p). Next, whenever v*B ^ 0, one can take K G L(— I p ) 
"sufficiently large", so that v*Wv > 0. Now, if v*B — 0, condition (b) guarantees 
that v*Wv > 0, so this part of the claim is established. 

(II) One can always write H — V*E V V for some non-singular V . The state- 
ment follows from applying the corresponding transformation of coordinates to the 
realization of F(s). 

(III) This follows from the proof of part (I) together with Proposition 14.21 noting 
that by assumption HA + A*H G P„, which is equivalent to AH^ 1 +H~ 1 A* G P„ . 
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(IV) See part (I) of Proposition I7.6| so the claim is established. □ 



Example 8.2. We here illustrate Proposition 18.11 by examining sets of rational 
functions QP(3,v, 1) associated with h(E„), with v = 0,1,2. Recall that, in 
most parts of this work, minimality of the realization is not assumed. 

(a a 
a a 6, with a, 61,62 real parameters. It 
b t -b 2 0/ 

realizes the function f (s) — ( b i~ b 2)( s ~ a ) 

First, for > a, b\ > 6|, f e GV min (3, 2, 1) C V, for a > 0, 6| > b\ , / € ^ mi „(3, 0, 1) 
and for a{b\ - b\) = 0, ^59(3,1,1). 

Next, if > a(&2 — b\) f is not generalized positive, but we show that one can 
always find a static output feedback of the form (|8.2p so that f c \ e QV(3, 1, 1). 
Indeed consider part II of the claim with £1,2- Condition 11(a) is satisfied. To 
check condition 11(b) note that the null-space of B* is given by vectors of the 
form Thus this condition now reads, 2a(6^ — b\) > as required. 

/ a + kb^ a, — kbib2 b]\ 

Indeed, the closed loop system matrix is L c \ —l a + kb 1 b 2 a-kb\ b 2 ) and for 

bi / ±62, it corresponds to QV m i n (3, 1, 1) functions whenever , a > k. For 
instance, for a = 1, 61 = 1, 62 = 1 and fc = — 1, one obtains L c i = L a from 
Example EUJ □ 

9. CONCLUDING REMARKS 

As it is often the case, the introduction of a novel concept opens the door to new 
research questions. Concerning the set QV{r, v, p) we here mention a sample of 
four problems. 

1. LMI techniques 

A comprehensive survey of the LMI approach to the PRL appeared in [29]. As 
already mentioned Proposition 17.61 suggests that large parts of LMI techniques can 
be extended to QV functions. For example, the use of LMI to render a closed loop 
function positive, through static output feedback, was addressed in literature, see 
e.g. [3HI33]. This can be extended in the spirit of Section [S] 

2. Lyapunov Order 

Recall that in Examples 15.11 and 17.41 we considered the system matrices L 1 and 
— L^ 1 and the respective C77^(3, 1, 1) rational functions ip 7 (s) — and 

V^( s ) = jfp Consider now L x := |(£ 7 + L^) =( 1 J and the associated 



rational function "4> x (s) = . In |16) a (partial) Lyapunov order was intro- 

duced in which L 7 < L x . The Lyapunov order was recently used in [42] . It 
is of interest to find an interpretation of this partial order in the framework of the 
rational functions ip~/( s ) and ipx(s)- 

3. Model order reduction 

One can exploit the convex structure of all system matrices associated with the 



^meaning that whenever L 1 g h(H), for some non-singular Hermitian H, it implies that for 
the same H, L x e L(H). 
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set QV(r,v,p) to try to introduce a scheme of model order reduction of uncertain 
systems in the spirit of I 1 Section 5]. 

4. Realization of Even and Odd QV functions. 

Recall that in the scalar case Odd functions map iR to itself while Even QV 
functions map iM. to R + . Both sets were addressed in [5] in the framework of 
rational functions. One can study properties of all Even and Odd functions within 
a prescribed set GV(r, u, p). 

For example, if one considers (for simplicity only real) realizations of functions in 
QV(3, 1, 1), the Odd and the Even cases can be parameterized by, 

/O a &i \ / — ai a 2 b\ ai+a 2 >0 

Lodd = [ a b 2 L evcn — a 2 ai b \ b#0 

\6i -6 2 J \ b -b d J d>0. 

_ (b\-bl)s i , s _ 6^(01+02) , j 



1podd{s) = VWen(s) 



It is interesting to note that in the framework of the associated Lyapunov equation 
(|1.5p . in the Odd case Q = 0, while in the iJwera case Q is diagonal. 
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